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SUMMARY 

Charts  are  derived  for  the  determination  of  the  power  spectrum  parameters  A and  NJ  i’-Jt  ^ 
for  vertical  acceleration,  at  an  aircraft's  centre  of  gravity,  due  to  atmospheric  turbulence. 

The  method  takes  account  of  spanwise  variations  in  gust  loading,  and  so  overcomes  the 
paradox  of  an  infinite  Nt)  which  is  found  with  a one-dimensional  gust  model. 
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Oylow 

Coefficient 

Reference  time 

Aspect  ratio 

Coefficient 

Lift  coefficient 

Theodorsen  function 

Wing  chord 

Coefficient 

Coefficient 

Factor  to  allow  for  transverse  variation  in  wing  loading 

Operator  denoting  Fourier  transform 

Frequency  cjs 

Interpolation  fraction 

Wing  loading  at  transverse  section  y 

Transfer  function 

Transfer  function  for  1-D  gust  field 
Transfer  function  for  2-D  gust  field 
Hankel  function 
Bessel  function 

Non-dimensional  measure  of  A 

Non-dimensional  measure  of  Nq 

Integral  scale  of  turbulence 

Rate  of  level  crossings  of  mean  value 

Power  spectral  density  of  w 

Power  spectral  density  of  y 

Wing  area 

Sears  function 

Distance  travelled  in  time  t 

Time 

Longitudinal  velocity  component  of  aircraft 

Aircraft  weight 

.!F(av) 

Vertical  component  of  gust  velocity 
Bessel  function 
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General  response  of  aircraft,  and  in  particular  y is  often  taken  to  be  the  vertical 
acceleration  at  the  aircraft  centre  of  gravity 


Transverse  co-ordinate 

3F{z) 

Elevation  of  aircraft 

Angle  of  attack 

Coefficient 

Coefficient 

Coefficient 

Coefficient 

Coefficient 

Aircraft  mass  parameter 
Air  density 

Distance  travelled  in  time  t 
Standard  deviation  of  w 
Standard  deviation  of  y 
Time 

Wagner  function 
Kiissner  function 
Non-dimensional  frequency 
Cut-off  frequency 


Note:  A dot  denotes  differentiation  w.r.t.  time 


1.  BSTRODUCnON 


The  power  spectral  method  of  analysing  gust  loads  on  an  aircraft  (see  ref.  I)  requires  the 
evaluation  of  two  parameters,  A and  Ng.  If  jit)  is  the  output  of  some  transducer  (in  the  present 
case  we  consider  the  vertical  acceleration  measured  at  the  aircraft’s  centre  of  gravity),  w(t)  is 
the  vertical  gust  velocity,  PJif)  is  the  power  spectrum  of  the  gust  velocity  and  P^^f)  is  the 
power  spectrum  of  y(/),  we  define 
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where  Pw  (/)  is  assumed  to  follow  the  von  Karman 
expression : 

-n  V [l+(a/)^]“'* 


a = 1 -339  X 2w  X — 


where 


(1.2) 


(1.3) 
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U is  the  aircraft  forward  velocity  and 
L is  the  integral  scale  of  turbulence; 

asymptotically  Pvif)  behaves  as 

The  transfer  function  of  the  aircraft  is 

' H{f)  = 9^[y{t)\l^[w{t)\  (1.5) 

where  ^[. . .]  denotes  the  Fourier  transform  of  the  quantity  in  brackets.  We  have 

P^f)=\H(f)\^Puif)  (1.6) 


The  transfer  function  may  be  derived  by  obtaining  the  response  of  the  aircraft  to  a 
sinusoidal  gust  of  arbitrary  frequency.  For  the  simple  case  of  an  aircraft  with  one  degree  of 
freedom  (heave),  and  a one  dimensional  variation  of  gust  (gust  velocity  varies  sinusoidally  along 
the  flight  path  but  is  invariant  in  the  spanwise  or  vertical  direction)  the  transfer  function  is  well 
known.  It  is  given  (see  Appendix)  by 
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where  the  Sears  function,  S(ai),  is  a complex  expression  which  represents  the  unsteady  wake 
effect  caused  by  gust,  and  the  Theodorsen  function,  C(«>),  is  another  complex  expression  repre- 
senting the  unsteady  wake  effect  caused  by  a change  in  the  angle  of  attack.  (Note  that  a vertical 
velocity  of  the  aircraft  produces  an  aerodynamic  effect  equivalent  to  a change  in  angle  of  attack.) 
The  parameter  w is  the  reduced  frequency 
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measured  in  radians  per  semi-chord,  and  is  the  mass  parameter 


2(ff'/S) 

pcg(dCLlda) 


(1.9) 


The  transfer  function  (1.7),  in  combination  with,  say,  the  von  Karman  power  spectrum 
leads  to  an  infinite  value  for  A'q.  This  is  basically  because  of  the  limitations  of  the  one- 
dimensional gust  model,  as  the  high  frequencies  which  contribute  largely  to  the  numerator  of 
equation  1.2  are  not  realistically  modelled  by  a gust  which  is  invariant  across  the  span.  In  the 
next  section  we  consider  a gust  structure  which  varies  sinusoidally  in  the  spanwise  direction  as 
well  as  the  axial  direction.  At  high  wave  numbers  (isotropic  turbulence  is  simulated  by  making 
the  wave  numbers  in  the  axial  and  spanwise  directions  equal)  the  upgusts  and  downgusts  on  the 
wing  cancel  and  a realistic  value  is  obtained  for  Nq. 


2.  THE  EFFECT  OF  SPANWISE  VARIATIONS  IN  THE  GUST  STRUCTURE 
Consider  a gust  which  varies  sinusoidally  across  the  span  according  to  the  equation 
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(This  equation  involves  the  simplification  that  the  gust  is  symmetrical  across  the  span,  but  this 
simplification  is  without  loss  of  generality  as  rolling  motions  caused  by  an  asymmetry  would  not 
affect  the  vertical  acceleration  at  the  centre  of  gravity.) 

Using  the  “strip  theory”  assumption  that  each  chordwise  element  behaves  as  a two- 
dimensional  aerofoil,  we  obtain  the  total  lift  on  the  wing  as  the  lift  corresponding  to  the  updraft 
at  the  centre-section,  multiplied  by  the  factor 

/'span  2 

F(a))=  g(>’)cos  ^j^dy  (2.2) 

span  U 

J -span/2 


where  g(y)  is  the  wing  loading  at  the  section  y.  The  precise  form  of  the  function  g(y)  is  only  of 
secondary  importance,  and  we  have  made  the  simple  assumption  that  the  wing  loading  at  any 
section  is  proportional  to  the  chord  at  that  section.  We  have  also  assumed  that  the  wing  plan 
form  varies  from  a delta  wing  at  aspect  ratios  of  two  or  less,  through  trapezoidal  to  rectangular 
at  aspect  ratios  of  eight  or  greater.  The  formula  for  the  chord  of  a trapezoidal  wing  is 


chord  = mean  chord 


(2.3) 


where  y is  the  spanwise  co-ordinate  measured  from  the  centre  of  the  fuselage,  and  y varies  from 
0 at  an  aspect  ratio  of  8 or  larger,  to  4 (corresponding  to  a delta  wing)  at  an  aspect  ratio  of  2 or 
less,  according  to  the  formula 
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Accordingly  the  factor  in  equation  2.2  becomes 


(2.4) 


F(<u)  = 


' I 

span  L \4  span/ 

J -span  2 


1 sin 
i>/R 


2nfy 

U 


dy 


>|’sin  (u»/R/2)T 

(cuAl/2)  J 


(2.5) 

(2.6) 


2 


and  so  for  trapezoidal  or  rectangular  wings  (y  4),  lf{<o)l^  behaves  asymptotically  as  to*^, 
whilst  for  delta  wings  (y  = 4),  behaves  asymptotically  as  w-*. 

Still  using  the  “strip  theory”,  we  can  now  write  the  transfer  function  for  two-dimensional 
gusts  as 

(2.7) 

= |//o(a))|^  . 1S(«.)|^  |f(a.)i^  (2.8) 

Figures  3 and  4 show  graphs  of  the  transfer  functions  |//o(/)i^,  '//i(/)F,  and  i//2(/)i^ 
together  with  input  and  output  spectra.  Figure  4 differs  from  Figure  3 only  in  the  formulae 
used  to  evaluate  the  effects  of  the  non-steady  aerodynamics.  Figure  3 is  based  on  the  exact 
equations  A.20  and  A.21  (see  Appendix),  whereas  Figure  4 is  based  on  the  fitted  curves,  A. 20a 
and  A.21a,  which,  it  may  be  seen,  depart  from  the  former  equations  at  non-dimensional  frequen- 
cies above  10. 


3.  THE  EVALUATION  OF  A AND  No 

The  integrals  which  arise  in  a numerical  evaluation  of  A and  No  are 
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(3.1) 
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and  the  first  of  these  will  converge  to  a limit  with  increasing /c  because  the  integrand  converges 
much  faster  than  /"‘.  However  the  integrand  in  the  second  integral  will  only  give  rise  to  a 
convergent  expression  if  the  transfer  function  corresponding  to  a two-dimensional  gust  structure 
is  used;  the  simple  one-dimensional  gust  model  gives  rise  to  a non-convergent  (infinite)  integral. 
The  fact  has  sometimes  led  to  an  implicit  assumption  that  No  is  theoretically  infinite,  and  that 
some  arbitrary  procedure  is  necessary  to  obtair  a realistic  finite  value  of  No-  Recently 
Houbolt  (ref.  2)  has  given  a method  for  evaluating  No  based  on  a two-dimensional  gust 
structure  model,  but  earlier  work  (Houbolt  (refs.  3,  4))  in  which  a design  manual  for  vertical 
gusts  had  been  given,  was  based  on  a one-dimensional  gust  model.  There,  it  was  recommended 
that  a cut-off  frequency, /f,  be  chosen  such  that  equation  3.1  was  within  a small  fraction  (2%) 
of  its  limit  value,  and  that  the  same  cut-off  frequency  be  applied  to  equation  3.2.  It  so  happens 
that  this  procedure  gives  approximately  the  right  answer  for  the  most  common  cases  with  an 
aspect  ratio  around  4 to  8 and  the  mass  parameter,  n,  in  the  range  10  to  50.  Houbolt  remarks 
that  the  No  value  derived  from  the  model  should  be  increased  by  a factor  of  the  order  of  2 to 
allow  for  flexibility  effects,  and  similarly  the  A value  should  be  increased  by  about  10"„.  Houbolt 
states  “this  factor  may  be  adjusted  upward  or  downward  if  some  additional  insight  relative  to 
the  aeroplane  response  characteristics  is  on  hand".  The  method  proposed  here  takes  account 
of  aspect  ratio  effects  in  a definite  way  which  would  otherwise  be  merely  left  to  the  engineer's 
judgement. 

The  transfer  function  of  the  aircraft  vertical  acceleration  from  two  dimensional  vertical 
gusts  is  a function  of  the  aspect  ratio,  /R,  and  the  non-dimensional  frequency,  w,  where 
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However,  the  von  Karman  power  spectrum  equation  is  a function  of  a different  form  of  non- 
dimensional  frequency,  a/,  where 
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and  so  the  solutions  to  equations  1.1  and  1.2  will  have,  as  parameter,  2Ljc.  Using  equations  2.7, 
2.6  and  A. 29  we  may  conveniently  obtain  the  solutions  in  the  non-dimensional  form 
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The  convergence  of  KJfi  and  towards  a limit  as  the  cut-off  non-dimensional  frequency, 
(Of,  is  raised  as  shown  in  Figure  5.  It  may  be  seen  that  A (or  KJix)  attains  its  limit  (or  nearly  attains 
i'  ')  at  a slightly  lower  cut-off  frequency  than  does  Nq  (or  kg). 


"USSION 

-g  equations  3.7  and  3.8  charts  have  been  derived  for  the  determination  of  A (see  Fig. 
ind  Nq  (see  Fig.  7).  Houbolt’s  chart  for  A is,  as  would  be  expected,  very  similar  to  the  charts 
presented  here  for  the  low  aspect  ratio  case.  There  is  more  variation  in  the  charts  for  A'q,  and 
Figure  8 shows  how  both  A and  Nq  are  affected  by  variations  in  aspect  ratio.  A is  little  affected 
except  at  large  aspect  ratio  and  low  values  of  fi.  However  Nq  is  strongly  affected  by  aspect  ratio, 
and  this  variation  seems  to  be  nearly  independent  of  n and  2L/c  except  at  quite  large  aspect  ratio. 

The  method  proposed  here  is  most  likely  to  be  in  error  at  low  aspect  ratios  where  the  "strip 
theory”  is  at  best  uncertain,  but  the  following  example  (Table  1)  shows  that  even  for  an  aspect 
ratio  as  low  as  2 the  method  produces  reasonable  agreement  with  observation. 
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TABLE  1 


Basic  Data 

Aircraft 

Mirage  IIIO 

Altitude 

600  m 

Speed 

232  m/s 

Aspect  ratio 

1 -97 

Aerodynamic  mean  chord 

5-25  m 

32 

2L/C 

88 

0 022 

A 

01  ® 
m/sec 

Determination  of  No  by  various  Methods 

(a)  One-dimensional  gust  model  (as  modified  by 

ko  = 0-073 
A'o  = 1 Hz 

(h)  Two-dimensional  gust  model  (Fig.  7) 


(r)  Experimental- 
(d)  Experimental- 


Rigid  aircraft 
-Flexible  aircraft 


ko  = 0-13 
A^o  = 1 -8  Hz 
Ao  = 2 Hz 
No  = 3-4  Hz 


For  this  example,  the  experimental  determination  of  Nq  has  been  described  by  Sherman 
(ref.  6).  In  brief  the  power  spectrum  was  calculated  out  to  a Nyquist  frequency  of  30  Hz,  using 
a 100  second  record  of  vertical  acceleration  during  fairly  severe  turbulence.  (Part  of  flight  2-135 
as  described  by  Rider  et  al.  (ref.  7).)  The  principal  structural  flexibility  effect  was  due  to  a funda- 
mental wing  bending  mode  at  about  II  Hz.  The  Nq  value  was  determined  from  equation  1.2 
modified  to  have  a finite  cut-off  frequency.  For  the  flexible  aircraft  A'o  was  determined  by 
integrating  right  up  to  30  Hz,  whilst  an  estimate  of  A^o  for  ^ rigid  aircraft  was  obtained  by 
reducing  the  cut-off  frequency  to  9 Hz,  which  eliminated  most  of  the  fundamental  structural 
resonance  mode.  The  rigid  (one-dimensional  gust)  model  of  Houbolt  estimates  an  Nq  of  only 
1 Hz,  whereas  the  rigid  (two-dimensional  gust)  model  on  which  Figure  7 is  based  estimates 
an  A^o  of  1 - 8 Hz,  which  compares  well  with  the  experimental  determination  of  2 Hz  for  a rigid 
aircraft.  The  effect  of  structural  flexibility  increased  Nq  to  3-4  Hz,  which  agrees  w'ell  with 
Houbolt’s  suggested  100"o  increase  to  allow  for  structural  flexibility. 

This  one  set  of  results  does  not  of  itself  validate  the  strip  theory  at  low  aspect  ratios,  but  it 
does  provide  support  for  the  method  proposed  herein  for  estimating  Nq- 


5.  CONCLUSION 

The  charts  for  determination  of  A and  Nq  presented  here  appear  to  give  more  realistic  values 
of  Nq  than  some  of  the  other  charts  presently  in  use. 
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APPENDIX 

The  Transfer  Function  for  1-D  Gusts 

Consider  a rigid  aircraft,  constrained  against  rotation  or  lateral  motion,  flying  at  a constant 
speed,  U,  and  altitude,  z,  with  a lift  force  L = \pU^SCl.  If  the  aircraft  now  encounters  a vertical 
gust  velocity,  Wg,  and  rises  with  a vertical  velocity,  i,  the  pseudo-steady  (i.e.  if  dynamic  effects 
were  to  be  neglected)  increment  in  lift  would  be 

= Lift  due  to  gust— Lift  due  to  vertical  motion 

The  two  components  in  this  equation  are  both  modified  by  the  dynamic  response  of  the  system, 
but  each  in  a different  way  because  the  change  in  incidence  due  to  gust  is  felt  first  at  the  leading 
edge  of  the  wing,  and  progresses  across  the  wing  at  speed  C,  whereas  the  change  in  incidence 
due  to  aircraft  vertical  motion  affects  the  whole  wing  simultaneously. 

When  an  aerofoil  changes  its  incidence  by  an  amount  Sa  at  time  / = 0,  th,.  extra  lift  SCz, 
builds  up  gradually  to  its  final  value,  SCloo.  We  denote  by  s the  (non-dimer  ional)  distance 
travelled  since  the  time  / = 0,  in  units  of  the  semi-chord,  so  that 


We  have 


(A.2) 


SC/,  = SC/,00  . <I>(s)  (A. 3) 

where  <I)(s)  is  known  as  the  Wagner  function  in  the  case  of  the  two-dimensional  wing,  and, 
according  to  Jones  (ref  8)  may  be  approximated  in  the  three-dimensional  case  by  the  formula 


j < 0:  <I>(5)  = 0 

s > 0:  <I>(s)  = 1 —D  . exp  (— Ss)— £ . exp  (— es) 


(A.4) 


where  D.  8,  E and  e are  functions  of  aspect  ratio. 

Similarly  when  an  aerofoil  penetrates  a unit  step  gust  at  / = 0,  the  lift  builds  up  as 


SCl  = SC/.oo'I'(5) 


(A.5) 


where  is  known  as  the  Kiissner  function  in  the  case  of  the  two-dimensional  wing  and  may 
be  approximated  in  the  three-dimensional  (finite  span)  case  by  the  formula 


j < 0:  T(s)  = 0 

s > 0:  TXr)  = I —A  . exp  (— a.r)  — 6 . exp  (— f3.s) 


(A.6) 


where  A,  a,  B and  ^ are  functions  of  aspect  ratio,  as  in  the  following  table  (see  Zbrozek  (ref  9) 
and  Bisplinghoff  et  at.  (ref  10)); 


Aspect 

Ratio 

Kiissner  function  'F(i) 

Wagner  function  0(5) 

A 

a 

B 

P 

D 

S 

E 

€ 

00 

lEl 

013 

mgm 

1-0 

0-335 

6 

0-294 

0-965 

0-361 

0-381 

3 

m 

0-558 

0-227 

3-20 

0-283 

0-54 

0 

For  calculation  purposes  we  need  an  interpolation  formula  to  give  the  value  of  <I>  or  T at 
other  aspect  ratios  than  the  ones  given  in  the  table.  If  ^411  and  /R2  are  two  of  the  aspect  ratios 
for  which  the  <I>  and  T functions  are  known,  and  ^41  is  some  intermediate  aspect  ratio  with 
/R1  < yR  < /R2,  then  the  interpolation  formula  used  here  is 


Vl/zll  - -v/1/^2 

y/\im  - 


(A.7) 


To  extrapolate  the  values  of  <1>  and  T for  aspect  ratios  below  3 (it  should  be  noted  that  a strip 
theory  such  as  we  will  use  later  is  not  an  accurate  approximation  at  such  low  aspect  ratios)  we 
have  extended  the  above  table  by  assuming  that  for  aspect  ratios  of  0 - 1 or  less  the  values  of  A, 
B,  D and  E attain  their  limiting  values  of  zero. 

The  theory  on  which  the  expressions  for  and  TCj)  have  been  derived  is  a linear  one,  so 
the  effect  of  complex  motions  can  be  built  up  by  a process  of  superposition.  The  result,  as  we 
now  show,  is  a convolution  (also  called  Duhamel  or  superposition)  integral.  If  the  vertical  air- 
craft velocity  changes  from  i to  z-f  Si  whilst  the  aircraft  position  changes  from  i = <t  to 
s = a+Sa,  the  increment  in  lift  at  some  later  position,  s (see  Fig.  2),  is 

SCt  = 8CLao‘I>(i-<r)  = 

dot  U 


*0  the  total  increment  in  lift  since  time  / = 0 is 

/*» 


ACl  = 


dCL  dz 
doL  da  U 


and  with 


we  obtain 


Tfi 


dz( 

dr 


<t>(s—a)da 


Tfl 


and  similarly  for  changes  Su’^  in  the  vertical  gust  velocity  we  obtain 


~ U da 


(JT  L C 


d 0 


(A.8) 


(A.9) 


(A.  10) 


(A.ll) 
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The  convolution  integrals  (equation  A.  10  and  A.  II)  can  conveniently  be  solved  by  Fourier 
transformation,  and  for  this  purpose  we  will  need  the  Fourier  transforms  (possibly  requiring  to 
be  understood  in  the  sense  of  generalised  functions)  of  the  functions  <l>  and  *F.  To  this  end,  con- 
sider the  case  of  the  harmonic  oscillations 


z = Zq  exp  (InifT) 

(A.  12) 

or 

= Wo  exp  (Inifr) 

(A.  13) 

The  solutions  to  these  two  cases  have  been  given  by  Sears  (ref.  11)  as 

ACt  = 2w[C(<o)4- Jioj]  exp  (2wift) 

(A.  14) 

and 

ACi  = 2nS(<o)  exp  (2nift) 

(A.15) 

respective,  where 

S((o)  = lJg(a))—iJ^(w)]C^w)-\-iJ^{w) 

(A.  16) 

//,'^'(a,)  + ,77o'^V) 

(A.  17) 

(A.  18) 

in  which  S(<o)  has  become  known  as  the  Sears  function,  and  C(<u)  is  the  Theodorsen  function. 

is  the  Hankel  function  of  the  second  kind,  and  to  is  the  non-dimensional  frequencv 
based  on  the  wing  semi-chord. 

2nf(cl2) 

ti>  = - 

U 

(A.  19) 

Sears  considered  a Joukowski  aerofoil  for  which  dCilda  = 2w.  With  this  substitution  we  may 
equate  ACt  from  equation  10  or  II  with  the  value  of  ACt  in  equation  14  or  15.  Then  using 
equations  12  and  13  and  taking  the  (formal)  Fourier  transform  we  obtain 


C(to)  4 J/to 

'2nif 

(A.20) 

S(<o) 

~ imf 

(A.21) 

The  first  of  these  equations  differs  a little  from  the  result  usually  quoted,  in  the  inclusion  of  the 
term  Jf<u,  which  includes  the  effect  of  virtual  mass  (the  so-called  non-circulating  lift]. 

As  an  alternative  to  equations  A. 20  and  A. 2 1 we  may  use  the  numerical  approximations  in 
equations  A.4  and  A. 6 to  obtain 


#•[<!)] 


1 

J 0^ 


txp  {—lirift)  dt 


D _ £ 

I + J “ • +”^j 

/a»  i(o ) 


(A.20fl) 
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(A.21a) 


However,  using  these  equations  it  will  be  necessary  to  account  separately  for  the  effect  of  virtual 


Equation  A.l  can  now  be  modified  to  take  account  of  dynamic  effects: 


AZ,  = Lift  due  to  gusts— Lift  due  to  a/c  vertical  motion 


aoL  jj  ar 

\^J  0 


and  the  vertical  equation  of  motion  may  be  written 


W 

^L  = —z 
g 


(A.23) 


which  with 


pcgdCildoL 


(A.24) 


becomes 


. , ''diw.rjt/,  ,1  f „p[/, 

Jo  *”  “*  Jo  *” 


(A.25) 


Denoting  by  Z and  Wg  the  Fourier  transforms  of  z and  Wg  respectively,  this  equation  becomes, 
under  Fourier  transformation 


if)^+i2-nifY3^m\Z  = InifWg^^] 


(A.26) 


and  since  the  Fourier  transform  of  the  aircraft  acceleration  is 


^[z]  = (ImffZ 


the  transfer  function  of  the  aircraft  becomes,  for  two-dimensional  gusts. 


H,{f)  = 


^(z)  ilni/WlT] 


Wg  f^  + 

U 


(A.27) 


-{rh 


(A.28) 


2tiwi+i27rif)^[<^] 


On  substituting  the  expressions  A. 20  and  A. 21  we  obtain 


^^y(2«>)^|5(a>)|' 


(A.29) 


|(2/i+i)/<u+C(<u)|^ 


= |//o(/)l'|5(o>) 


(A.30) 


J 


Hoif)  = 


(A.31) 


(2/i + \)i<x> + C(«o) 


is  the  transfer  function  obtained  by  neglecting  the  effect  of  unsteady  aerodynamics  of  gust 
penetration  (but  not  of  the  less  important  effect  of  unsteady  aerodynamics  of  the  vertical  motion 
of  the  aircraft). 

It  may  be  seen  that  the  effect  of  the  virtual  mass  is,  as  would  be  expected,  the  same  as  a slight 
increase  in  the  mass  parameter,  m,  equivalent  to  increasing  the  aircraft  mass  by  the  mass  of  a 
cylinder  of  air  of  diameter  equal  to  the  wing  chord. 

The  function  C(iu)  behaves  asymptotically  as  a non-zero  constant,  |5{<o)p  behaves  as 
l/(l  + 2n-<u)  (the  Liepman  approximation)  and  so  |//o(/)|^  behaves  asymptotically  as  w-'. 


I 
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FIG.  1 CHANGE  IN  ANGLE  OF  ATTACK  DUE  TO  GUST  VELOCITY 
AND  AIRCRAFT  VERTICAL  MOTION 
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FIG.  3 INPUT  SPECTRUM,  TRANSFER  FUNCTION  AND  OUTPUT 
(RESPONSE)  SPECTRUM  (USING  EQUATIONS  A.20,  A.21) 
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FIG.  4 INPUT  SPECTRUM,  TRANSFER  FUNCTION  AND  OUTPUT 
(RESPONSE)  SPECTRUM  (USING  APPROXIMATE  EQUATIONS  A.20a,  A.21a) 
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FIG.  5 CONVERGENCE  OF  THE  INTEGRALS  FOR  AND  ko  TOWARDS 
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FIG.  6 CHARTS  FOR  THE  DETERMINATION  OF  Nq  FOR  c.g.  VERTICAL  ACCELERATION 
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